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>• ' Abstract 

o : 

I In this lecture I review the theory of magnetized accretion discs driving jets, with 

00 ' a focus on Young Stellar Objects (YSOs). I first introduce observational and theoreti- 

, cal arguments in favor of the "disc wind" paradigm. There, accretion and ejection are 

interdependent, requiring therefore to revisit the standard picture of accretion discs. A 
simple magnetostatic approach is shown to already provide some insights of the basic 
. phenomena. The magnetohydro dynamic equations as well as all usual assumptions are 

\^ I then clearly listed. The relevant physical mechanisms of steady-state accretion and ejec- 

tion from Keplerian discs are explained in a model independent way. The results of 
self-similar calculations are shown and critically discussed, for both cold and warm jet 
i configurations. I finally provide observational predictions and the physical conditions re- 



quired in YSOs discs. The necessity of introducing a magnetospheric interaction between 



p i' the disc and the protostar is briefiy discussed. 

o 

^ : 1 Introduction 

Collimated ejection of matter is widely observed in several astrophysical objects: inside 
^ I our own galaxy from all young forming stars and some X-ray binaries, but also from the 

' core of active galaxies. All these objects share the following properties: jets are almost 

^ cylindrical in shape; the presence of jets is correlated with an underlying accretion disc 

surrounding the central mass; the total jet power is a sizeable fraction of the accretion 
power. 



1.1 Why should jets be magnetized ? 

There are basically two kinds of jet observations: spectra (blue shifted emission lines) and 
images (in the same lines for YSOs, or in radio continuum for compact objects). Most 
of these images show jets that are extremely well collimated, with an opening angle of 
only some degrees. On the other hand, the derived physical conditions show that jets 
are highly supersonic. Indeed, emission lines require a temperature of order T ~ 10^ K, 
hence a sound speed Cg ^ 10 km/s while the typical jet velocity is Vj ~ 300 km/s. The 
opening angle of a ballistic hydrodynamic fiow being simply tan^ = Cg/vj, this provides 
~ 5° for YSOs, nicely compatible with observations. Thus, jets could well be ballistic, 
showing up an inertial confinement. But the fundamental question remains: how does a 
physical system produce an unidirectional supersonic Row ? Naively, this implies 
that confinement must be closely related to the acceleration process. 
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This has soon been recognized in the AGN community (where jets were being observed 
for quite a long time, see Bridle & Perley 1984), leading Blandford & Rees to propose in 
1974 the "twin-exaust" model. In this model, the central source is emitting a spherically 
symmetric jet, which is confined and redirected into two bipolar jets by the external 
pressure gradient. Indeed, the rotation of the galaxy would probably produce a disc-like 
anisotropic distribution of matter around its center. Thus, in principle, the ejected plasma 
could be focused towards the axis of rotation of the galaxy (where there is less matter 
there) and thereby accelerated like in a De Laval nozzle. Such an idea was afterwards 
applied by Canto (1980) and Konigl (1982) for YSOs. However, this model had severe 
theoretical drawbacks, related to Kelvin-Helmholtz instabilities quickly destroying the 
nozzle. Besides, it docs not explain the very origin of the spherical flow. 

But the strongest argument comes from recent observations, like those of IIII30 (see 
eg. Ray et al. 1996). We know now that jets are already clearly coUimated close to the 
central star, with no evidence of any relevant outer pressure. This implies that jets must 
be self-collimated. In my opinion, such an observation also rules out the proposition that 
jets are collimated by an outer poloidal magnetic field (Spruit et al. 1997). 

The only model capable of accelerating plasma along with a self-confinement relies 
on the action of a large scale magnetic field carried along by the jet. In fact, Lovelace 
and Blandford proposed independently in 1976 that if a large scale magnetic field would 
thread an accretion disc, then it could extract energy and accelerate particules (electron 
positron pairs in their models). Then, Chan Sz Henriksen (1980) showed, using a simplified 
configuration, that such a field could indeed maintain a plasma flow collimated. But it 
was Blandford & Payne who, in 1982, produced the flrst full calculation of the interplay 
between a plasma flow (made of electrons and protons) and the magnetic field, showing 
both acceleration and self-collimation. 



1.2 What are the jets driving sources ? 

To make a long story short, there are three different situations potentially capable of 
driving magnetized jets from young forming stars^. 

• the protostar alone: these purely stellar winds extract their energy from the 
protostar itself (eg. Mestel 1968, Hartmann & McGregor 1982, Sauty et al. 1999). 

• the accretion disc alone: "disc winds" are produced from a large radial extension 
in the disc, thanks to the presence of a large scale magnetic field (eg. Blandford & 
Payne 1982, Pudritz &; Norman 1983). They are fed with both matter and energy 
provided by the accretion process alone. 

• the interaction zone between the disc and the protostar: these "X-winds" 
are produced in a tiny region around the magnetopause between the disc and the 
protostar (eg. Shu et al. 1994, Lovelace ct al. 1999, Ferreira et al. 2000). 

Purely stellar wind models are less favoured because observed jets carry far too much 
momentum. In order to reproduce a YSO jet, a protostar should be either much more 
luminous or rotating faster than observations show (DeCampli 1981, Konigl 1986). This 
leaves us with either disc- winds or X-winds^. Prom the observational point of view, it is 
very difficult to discriminate between these two models. In a nice review, Livio (1997) 

^An alternative model is based on some circulation of matter during the early infall stages (see Lery et al. 
1999 and references therein). However, by construction, such a model is only valid for Class sources and 
cannot be used to explain jets from T-Tauri stars. See also Contopoulos & Sauty (2001). 

^Camenzind and collaborators proposed an "enhanced" version of stellar winds, related to an old idea of 
Uchida & Shibata (1984). In this picture, a magnetospheric interaction with the accretion disc is supposed to 
strongly modify both jet energetics and magnetic configuration, leading to enhanced ejection from the protostar 
(Camenzind 1990, Fendt et al. 1995, see Breitmoser & Camenzind 2000 and references therein). 
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gathered a number of arguments for the disc wind model. The main idea is to look for a 
model able to explain jets from quite a lot of different astrophysical contexts (YSOs, AGN, 
X-ray binaries). The only "universal" ingredient required is an accretion disc threaded 
by a large scale magnetic field. Such a paradigm naturally explains (qualitatively) all 
accretion-ejection correlations known and is consistent with every context: young stars 
(see Cabrit's contribution, this volume), microquasars (eg. Mirabel & Rodriguez 1999) 
and AGN (see eg. Serjeant et al. 1998, Cao & Jiang 1999, Jones et al. 2000). 

1.3 Where does this magnetic field come from ? 

Let's face it: we don't know. There are two extreme possibilities. The first one considers 
that the field has been advected by the infalling material, leading to a flux concentration 
in the inner disc regions. The second one relies on a local dynamo action in the disc. 
Most probably, the answer lies between these two extreme cases. 

If the interstellar magnetic field has been indeed advected, the crucial issue is the 
amount of field diffusion during the infall. Indeed, if we take the fiducial values n ~ 1 cm~^ 
and S ~ 4 ^G of dense clouds and use the law B oc n^/^ (Crutcher 1999), we get a 
magnetic field at 1 UA ranging from 10 to 10'^ G! One must then consider in a self- 
consistent way the dynamical infiuence of the magnetic field, along with matter energy 
equation and ionization state (see Lachieze's contribution). This is something extremely 
difficult and, as far as I know, no definite result has been yet obtained (see however Ciolek 
& Mouschovias 1995). 

How exactly disc dynamo works is also quite unclear. Dynamo theory remains intri- 
cate, relying on the properties of the turbulence triggered inside the disc. In our current 
picture of accretion discs, these are highly turbulent because of some instability, probably 
of magnetic origin (sec Tcrquem's contribution). Such a turbulence is believed to provide 
means of efficient transport inside the disc, namely anomalous viscosity, magnetic diffu- 
sivity and heat conductibility. Obviously, small scale (not larger than the disc vertical 
scale height), time dependent magnetic fields will then exist inside accretion discs. But 
we are interested in the mean flow (hence mean field) dynamics. So, in practice what does 
ejection require ? To produce two opposedly directed jets, there must be a large scale 
magnetic field in the disc, which is open and of one of the following topologies: 

Dipolar: the field threads the disc, with only a vertical component at the disc mid- 
plane, matter being forced to cross the field lines while accreting (eg. Blandford & Payne 
1982). 

Quadrupolar: field lines are nearly parallel to matter, entering the disc in its plane 
and leaving it at its surfaces, with only a non zero azimuthal component at the disc 
midplane (eg. Lovelace et al. 1987). 

Most jet models and numerical simulations assume a dipolar magnetic configuration, 
with no justification. In fact, it turns out from the analysis of disc physics that only 
the dipolar configuration is suitable for launching jets from Keplerian accretion discs (see 
appendix A in Ferreira 1997) . This has been recently confirmed using a dynamo-generated 
magnetic fields (Rekowski et al. 2000). This is due to a change of sign of the a effect across 
the disc midplane, as observed in numerical simulations of MHD turbulence in the shearing 
box approximation (Brandenburg Sz Donner 1997). But a realistic situation requires 
to treat both turbulence and the backreaction of the magnetic field in a self-consistent 
way. Indeed, jet production is a means of ffux leakage, hence of possible dynamo self- 
regulation (Yoshizawa &; Yokoi 1993). Anyway, this severe issue of dynamo and turbulence 
lead theorists to simply assume the existence of a large scale magnetic field. Its' value 
and distribution are then either imposed or obtained as conditions for stationarity in 
Magnetized Accretion-Ejection Structures (MAES). 
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1.4 The MAES paradigm 



A large scale (mean) magnetic field of bipolar topology is assumed to thread an accretion 
disc, allowing ejected plasma to flow along open field lines. This field extracts both angular 
momentum and energy from the underlying disc and transfers them back to the ejected 
plasma. There have been numerous studies of magnetized jets (e.g. Blandford & Payne 
1982, Heyvaerts & Norman 1989, Pelletier & Pudritz 1992, Contopoulos & Lovelace 1994, 
Rosso & Pelletier 1994, Lery et al. 1999 to cite only a few), but they all suppose that the 
underlying disc would support the jets. In all these works, the disc itself was treated as a 
boundary condition, people usually assuming a standard viscous accretion disc (Sliakura 
& Sunyaev 1973). However, if jets are to carry away the disc angular momentum, they 
strongly influence the disc dynamics. 

The investigation of MAES, where accretion and ejection are interdependent, requires 
a new theory of accretion discs. The relevant questions that must be addressed by any 
realistic model of stationary magnetized accretion-ejection structures are the following: 

(1) What are the relevant physical mechanisms inside the disc ? 

(2) What are the physical conditions allowing accretion and ejection ? 

(3) Can we relate jet properties to those of the disc ? 

To answer these questions, one must take into account the full 2D problem (not 3D, 
thanks to axisymmetry) and not treat the disc as infinitely thin as in a standard disc 
theory. As a consequence, no toy-model has been able yet to catch the main features 
of these accretion-ejection structures. There, the disc accretion rate exhibits a radial 
variation as matter is being ejected. The link between accretion and ejection in a MAES 
can therefore be measured by the quantity ^ = '^^"^f^" called the "ejection index". This 
parameter (which can vary in the disc) measures a local ejection efficiency. In a standard 
accretion disc, ^ = everywhere leading to a constant accretion rate. A complete theory 
of MAES must provide the allowed values of ^ as a function of the disc properties. 



2 Theoretical framework of MAES 

2.1 A magnet ost at ic approach 
2.1.1 The Bcirlow experiment 

A wheel made of a conducting material is put between the two poles of an electromagnet. 
This device produces a magnetic field perpendicular to the disc, along its axis of rotation. 
The wheel brushes against mercury contained in a tank, thereby allowing to close an 
electric circuit (mercury, wire connecting the disc axis, the disc itself). In order to check 
if some current is flowing, we can put a small lamp. With a crank, we then provide a 
rotating motion to the disc and let it evolve. If the magnetic field is off, nothing happens: 
the lamp stays mute (no current) and the disc very gradually slows down. On the contrary, 
if the magnetic field is present, the lamp lights on (a current is flowing) and the disc stops 
very quickly ! When the disc is finally motionless, the lamp is also off. 

The explanation of this phenomenon lies in electromagnetic induction. The disc is 
made of a conducting material, meaning that charged particules are free to move inside 
it. When the disc starts to rotate, these particules drag the magnetic field along with them. 
The field lines become then twisted, showing a conversion of mechanical into magnetic 
energy. Once the disc has stopped, all the initial mechanical energy has been converted 
(here, dissipated into heat along the whole electric circuit). The magnetic field was just 
a mediator between different energy reservoirs. 
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2.1.2 On the importance of currents 

A MAES is an astrophysical Barlow wheel. Indeed, the rotation of a conducting material 
in a static magnetic field induces an electric (more precisely electromotive) field 

Em = uABc=l Q.rBzer (1) 

where Q is the matter angular velocity and e^- the radial unit vector (in cylindrical coor- 
dinates). This field produces an electromotive force e, ie a difference of voltage between 
the inner and outer disc radii, namely 



rr 
Jn 



Em-dr- ^irfBi (2) 



(where the dominant contribution is assumed to arise at the inner radius). In steady 
state, two electric circuits can develop (above and below the disc midplane) with a current 
I± = e/Z± where Z± is their impedance. There is therefore an available electric power 

P = e(/_+J+) = i^ (3) 

where Zeq is the equivalent impedance. Since a current J = /_ + J+ is flowing inside the 

disc, it becomes affected by the magnetic field through the Laplace force F = J^" Idr AB, 
which provides a torque acting against the disc rotation. This result is consistent with 
Lenz's law, which states that induced currents work against the cause that gave birth to 
them. 

Now, the existence of such a torque allows matter to fall towards the central engine, 
hence to be accreted. Energy conservation implies that the mechanical power liberated 
through accretion, namely 

P (4) 

where Ma is the accretion rate, must be equal to the electric power P. This imposes an 
impedance matching 

2rfBf [gM 

where fio is the vacuum permeability. The rhs estimate used a magnetic field close to 
equipartition with the thermal energy inside the disc (see Sect 5.2.2). The disc accretion 
rate depends therefore on the global electric circuit. 

Electric currents I± flow along the axis, more precisely are distributed inside the jets. 
But if we model such current as being carried by an electric wire, we can estimate the 
generated magnetic field, namely -B<^(r) = — /Xo-^/27rr. Such a field is negative, consistent 
with the shear given by the disc rotation. Because of this field, there is a Laplace force 
towards the axis, ie the famous "hoop stress" which maintains the jets coUimated. 

To summarize: 

1. Accretion: a torque due to the magnetic field extracts the disc angular momen- 
tum. Such a torque is related to the establishment of two parallel electric circuits. 
Their equivalent impedance is directly linked to the accretion rate affordable by the 
structure. 

2. Global energy budget: a fraction of the accretion power is dissipated into Joule 
heating (disc and jets heating, allowing radiation) , another is converted into kinetic 
power carried by the jets. 
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3. Jet acceleration: possible due to the conversion of electric into kinetic power. 
Note that if the impedances and Z_ (which describe all dissipative effects) are 
different, then two bipolar but asymmetric jets can be produced. 

4. Jet collimation: jets do have naturally a self-coUimating force if they carry a non- 
vanishing current. Since the electric circuit must be closed, not all magnetic field 
lines embrace a non-zero current. This obviously implies that all field lines anchored 
in the disc cannot be self-coUimated (Okamoto 1999). 

Lots of physics can be understood within the framework of magnetostatics. However, 
the precise description of a rotating astrophysical disc, its interrelations with outflowing 
plasma as well as the calculation of the asymptotic current distribution inside the jets 
(necessary to understand collimation) quite evidently require a fluid description. 

2.2 Magnet ohydrodynamics 

Magnetohydrodynamics (MHD) describes the evolution of a coUisional ionized gas (a 
plasma) submitted to the action of an electromagnetic field. Because the source of the 
field lies in the motions of charged particules (currents), the field is intrinsically tied to 
and dependent of these motions. Therefore, because of this high non-linearity, MHD offers 
an incredible amount of behaviours. 

2.2.1 Prom a multicomponent to a single fluid description 

Circumstellar discs and their associated jets are made of dust and gas, which is composed 
of different chemical species of neutrals, ions and electrons. Wc should therefore make a 
multicomponent treatment. However, this is far too complex especially when energetics 
comes in. On the other hand, if all components are well coupled (through collisions) a 
single fluid treatment becomes appropriate. 

For each specie a, we define its numerical density n^, mass ttt-q,, electric charge Qa, 
velocity Va and pressure Pa- The equation of motion for each specie writes 

Pa^ = -^Pa - Pa^^G + ^ i^/3a + "^^^(^ + ^ B) (6) 
Ut ^ 

where D./Dt = d./dt + • V is the Lagrangean derivative, pa = Uama, is the 
gravitational potential of the central star and Fpa is the collisional force due to all other 
species We can define the "mean" flow very naturally as 

P = ^ Uania 
a 

PU = '^anaVa 
a 

P = ^naksT (7) 

a 

J = ^naQaVa 

a 

where p is the density, u the velocity, P the pressure, J the current density and ks the 
Boltzmann constant. A single fluid description becomes relevant whenever the plasma is 
enough collisional. In such a case, we can safely assume that all species share the same 
temperature T. Moreover, we assume that any drift between the mean flow and a specie 
a is negligible, namely \\va — u\\ ^ H""!!- Using Newton's second law {Y,a,pPai3 = 0) and 
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local electrical neutrality {J2a''^aQa = 0), we get the usual dynamical equations for one 
fluid 



^ + V-pu = (8) 
p— = -VP-pV^G + JAB (9) 

by summing over all species a. Even if the bulk of the flow is made of neutrals, they 
feel the magnetic force through collisions with ions (mainly) and electrons, J A B = 
(1 + X){Fin + Fen) where X = pi/pn is the density ratio of ions to neutrals. 

The evolution of the electromagnetic fleld is described by Maxwell's equations, namely 
(in vacuum) 

V-B = (10) 
1 dE 

VAB = ^oJ+^^ (11) 
V-E = ^ (12) 

VAE = -f (13) 

Faraday's law (Eq. 13) shows that the strength of the electric field varies like E/B Ljt ~ 

U where [7 is a typical plasma velocity. Now, using Ampere's equation (11), shows that 
the displacement current has an effect of order (C//c)^ only with respect to real currents 
(c is the speed of light). Thus, in a non-relativistic plasma, it can be neglected providing 
a current density 

J=— VAB (14) 



o 



directly related to the magnetic field. Under this approximation, electric charge conser- 
vation 

^ + V-J = (15) 

shows that no charge accumulation is allowed: the first term is also of order (U/cf'. 
Therefore V • J = 0, implying closed electric circuits. 

Energy conservation of the electromagnetic field writes 

dW ^ -> 

+ V ■ Smhd = -J E (16) 

where W = We + Wm = SoB'^/2 + B^/2po is the electromagnetic field energy density and 

Smhd = (17) 

po 

is the Poynting vector, carrying the field energy remaining after interaction with the 
plasma (the term J ■ E). Inside the non-relativistic framework, the energy density con- 
tained in the electric field is negligible with respect to Wm (of order {U/c)^). 

2.2.2 Generalized Ohm's law 

In order to close the above system of equations, we need to know the electric field E. 
Its' expression is obtained from the electrons momentum equation, consistently with the 
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single fluid approximation. Namely, we assume that electrons are so light that they react 
almost instantaneously to any force, i.e. 

= -VPe + Y^Fpe-ene{E + VeAB) (18) 

E + uAB = {iI-^^)AB-^ + y^^^{v3-Ve) (19) 

eue ^ eue 

Now, using (i) the expression of the Lorentz force JAB, (ii) the approximation J ~ 
erieivi — Ve) and (iii) neglecting the contribution due to the collisions between electrons 
and neutrals (with respect to those involving ions), we obtain the generalized Ohm's law 

^ - - JAB f pnV [J A B) A B VPe . X 

E + uAB = riJ + — ^ (20) 



eUe VP/ minriiVin eUe 

where r] = {ninenn^ne + "i-ie'T-j^ie)/(erae)^ is the electrical (normal) resistivity due to colli- 
sions, and niap the reduced mass. The first term on the rhs is the Ohm term, the second 
is the Hall effect, the third the ambipolar diffusion term and the fourth a source of elec- 
tric field due to any gradient of electronic pressure. Fortunately, all these terms become 
negligible with respect to uA B whenever the plasma is well coupled and ionized (p„ <S p) ■ 

2.2.3 Plasma energy equation 

This is the trickiest equation and can be written in several forms. We choose to express 
the internal energy equation, which can be written as 

,r- = r-A (21) 

where S is the plasma specific entropy, V all heating terms and A all cooling terms. Note 
that a transport term (eg. such as conduction) can cool the plasma somewhere and heat 
it elsewhere. The MHD heating rate due to the interaction between the plasma and the 
electromagnetic field writes 

- - - \2 JAB 
TMHD = J-{E + ilAB)c^riJ^+[i^] (22) 



The first term is the Joule eff'cct while the second is the heating due to ambipolar difi'usion. 
Note that, although dynamically negligible in discs, such an effect might be responsible 
for jets heating (Safier 1993, Garcia et al. submitted). 

2.3 Modelling a MAES 
2.3.1 Assumptions 

Our goal is to describe an accretion disc threaded by a large scale magnetic field of bipolar 
topology. In order to tackle this problem, we will make several symplifing assumptions: 

(i) Single-fluid MHD: matter is ionized enough and all species well coupled. As a 
consequence, we use the simple form of Ohm's law 

E + uAB = r]J (23) 

Taking the curl of this equation and using Faraday's law, we get the induction equation 
providing the time evolution of the magnetic field 

_ = VA(iiA^)- VA(:/VA5) (24) 
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where v = rj/fj,o is the magnetic difFusivity. The first term describes the effect of advection 
of the field by the flow while the second describes the effect of diffusion, matter being able 
to cross field lines thanks to difFusivity. 

(ii) Axisymmetry: using cylindrical coordinates (r, (p, z) no quantity depends on cp, 
the jet axis being the vertical axis. As a consequence = Q and all quantities can be 
decomposed into poloidal (the (r, z) plane) and toroidal components, eg. u = Up + QreJ 
and B = Bp + -B^eJ. A bipolar magnetic configuration can then be described with 

Bp = ^Va A , (25) 

where a(r, z) is an even function of z and with an odd toroidal field ^^(r, —z) = —B^{r, r). 
The flux function a is related to the toroidal component of the potential vector (a = rAfp) 
and a{r,z) = constant describes a surface of constant vertical magnetic flux $, 



I B-dS = 27ra(r, z) . (26) 

s 



Is 

The magnetic field distribution in the disc as well as the total amount of flux are unknown 
and must therefore be prescribed. 

(iii) Steady-state: all astrophysical jets display proper motions and/or emission 
nodules, showing that they are either prone to some instabihties or that ejection is an 
intermittent process. However, the time scales involved in all objects are always much 
larger than the dynamical time scale of the underlying accretion disc. Therefore, a steady 
state approach is appropriate^. 

(iv) Transport coefficients: if we use a normal (collisional) value for i/, we find that 
the ratio of advection to diffusion in (24), measured by the magnetic Reynolds number 

Tlm^^, (27) 

is always much larger than unity inside both the disc and the jets. This is the limit of ideal 

MHD where plasma and magnetic fields are "frozen in" . Jets are therefore described with 
ideal MHD [rj = 0). Within this framework, the stronger (initially the magnetic field) 
carries the weaker (ejected disc plasma) along with it. But inside the disc, gravitation is 
the dominant energy source and the plasma drags and winds up the field lines. Such a 
situation cannot be maintained for very long. Instabilities of different kind will certainly 
be triggered leading to some kind of saturated, turbulent, disc state (eg. tearing mode 
instabilities, or magneto-rotational instability, Balbus & Hawley 1991). 

In turbulent media, all transport effects are enhanced, leading to anomalous transport 
coefficients. These coefficients are the magnetic difFusivity Ufyi a,ncl resistivity Ty^, but also 
the viscosity (associated with the transport of momentum) and thermal conductivity kt 
(associated with the heat flux). Providing the expressions of these anomalous coefficients 
requires a theory of MHD turbulence inside accretion discs. Having no theory, we will 
use simple prescriptions, like the alpha prescription used by Shakura & Sunyacv (1973). 
In particular, because of the dominant Keplerian motion in discs, we allow for a possible 
anisotropy of the magnetic diffusivity. Namely, we use two different turbulent diffusivities: 

(related to the diffusion in the (r, z) plane) and (related to the diffusion in the (f) 
direction) . 

Since stationary discs must be turbulent, all fields (eg. velocity, magnetic field) must 

be understood now in some time average sense. 



^Note that this conclusion holds even in microquasars. In GRS 1915+105 the duration of an ejection event 
is around 10^ sec only, but the disc dynamical time scale is around 1 msec, ie. 10~^ times smaller. See however 
Tagger & Pellat (1999) for an alternative view on this topic. 
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(v) Non-relativistic MHD framework: in YSOs, matter remains always non- 
relativistic but there is something more about it. Indeed, magnetic field lines are anchored 
in an rotating accretion disc. Hence, if a field line of angular velocity Jl* opens a lot, there 
is a cylindrical distance such that its linear velocity reaches the speed of light, defining 
a "light cylinder" Rl = c/fi*. Now, if one imposes ideal MHD regime along the jet, 
ie. E = —u f\ B = — il^re^ /\ Bp, we see that the displacement current is no longer 
negligible at the light cylinder: propagation effects become relevant and one must take 
into account a local non-zero electric charge (provided by the Goldreich-Julian charge 
V ■ E = p^/eo)- As a consequence, the plasma feels an additional electric force p^E, even 
if the bulk velocity of the flow is non relativistic (eg. Breitmoser & Camenzind 2000)! 

But remember that this extra force and its corresponding "light cylinder" arose because 
of the assumption of ideal MHD. In fact, taking into account a local charge density 
probably imposes to also treat the non- ideal contributions (see Eq. 20). Nobody provided 
yet a self-consistent calculation. We just assume here that any charge accumulation would 
be quickly canceled (no dynamical relevance of the light cylinder). 



2.3.2 Set of MHD equations 

We use the following set of MHD equations: 
Mass conservation 

V • pn = (28) 

Momentum conservation 

pu-Vu = - VP - pV$G + JA -h V • T (29) 

where jipJ = V A S is the current density and T the turbulent stress tensorm which is 
related to the turbulent viscosity (Shakura & Sunyaev 1973). 
Ohm's law and toroidal field induction^ 

■nmJ(l> = Up A Bp (30) 
V-(%VrS<^) = V ■ -{B^up - Bpnr) (31) 



where r]m = p-oi^m and rj!^ = p-o^'m the anomalous resistivities. 
Perfect gas law 

P = p^^T (32) 
prup 

where nip is the proton mass and p a generalized "mean molecular weight" . 
Energy equation 

As seen previously, the exact energy equation (21) involves various physical mechanisms. 
Its explicit form is 

V-{UUp + Srad + Qturb) = -P^ ■ Up + r]m4 + Vm^p + VvlrVilf (33) 

— * 

where r)y = pp^, U is the internal specific energy, Sj-ad = —-^^Prad is the radiative energy 
flux (k is the Rosscland mean opacity of the plasma, Prad the radiation pressure) and qturb 
the unknown turbulent energy flux. This flux of energy arises from turbulent motions and 
provides both a local cooling Af^rb and heating Tfy^^b- Indeed, using a kinetic description 
and allowing for fluctuations in the plasma velocity and magnetic field, it is possible to 
show that all energetic effects associated with these fluctuations cannot be reduced to 
only anomalous Joule and viscous heating terms. Therefore, a consistent treatment of 



^Obtained from Eq. (24) and after some algebra (remember that E^ = Q). 
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turbulence imposes to take into account qturb (see early paper of Stiakura et al. 1978). 
But how to do it ? Moreover, the radiative flux depends on the local opacity k of the 
plasma, which varies both radially and vertically inside the disc. Besides, the expression 
of the radiative pressure Pj-ad is known only in optically thick media {Prad = while 
the disc can be already optically thin at the disc midplane. 

Thus, it seems that a realistic treatment of the energy equation is still out of range. 
As a first step to minimize its impact on the whole structure, we will use a polytropic 
equation 

P = Kp^ , (34) 

where the polytropic index F can be set to vary between 1 (isothermal case) and 7 = f 
(adiabatic case) for a monoatomic gas. Here K can be allowed to vary radially but remains 
constant along each field line. In section 5.3, we will turn our attention to thermal effects 
and use therefore a more appropriate form of the energy equation. 

Using the above set of MHD equations to describe astrophysical discs and jets would 
normally require consistency checks, at least a posteriori. If one fluid MHD seems justifled 
inside the inner regions of accretion discs, it may well not be anymore the case along the 
jet (huge fall in density). In particular non-ideal effects are most probably starting to play 
a role in Ohm's law (ambipolar diffusion or even Hall terms), allowing matter to slowly 
drift across field lines. This may have important dynamical effects dowstream. To clearly 
settle this question, the full thermodynamics of the plasma including its ionization state 
should be self-consistently solved along the jet. 



2.4 Critical points in stationary flows 

In the real world, everything is time dependent. Imagine matter expelled from the disc 
without the required energy: it will fall down, thereby modifying the conditions of ejection. 
A steady state is eventually reached after a time related to the nature of the waves 
travelling upstream and providing the disc information on what's going on further up. 

The adjustment of a MAES corresponds to the phenomenon of impedance matching. 
As we saw, this matching relates the accretion rate to the dissipative effects in the electric 
circuit. In practice, this means that the resolution of stationary fiows requires to take 
into account, in some way, all time-dependent feedback mechanisms. This is done by re- 
quiring that, once a steady-state is achieved, no information (ie. no waves) can propagate 
upstream, from infinity (in the z-direction) to the accretion disc. There is only one way 
to do it. Matter must flow faster than any wave, leaving the disc causally disconnected 
from its surroundings. 



2.4.1 The Parker wind 

Let us first look at the simple model of a spherically symmetric, isothermal, hydrodynamic 
fiow. Such a model was first proposed by Parker (1958) to explain the solar wind. In 
spherical coordinates, mass conservation and momentum conservation write 

^ + ^ + 2^ = (35) 
u p r 

du IdP GM 

u— H — H TT- = (36) 

dr p dr 

where gas pressure is P = pCg, the sound speed Cg being a constant. In such a simple 
system is hidden a singularity. Indeed, after differentiation one gets 



^ ^ ^ dr J. J.2 V ^ 
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{Cl-u^)—— = -2^ + ^ 38 
^ ^ ^ dr J. J.2 V / 

showing that the system is singular at r = r^, where u = Cs- In order to obtain a 
stationary solution, one must then impose the regularity condition which is 

GM 

In practice, fixing the distance of the sonic point imposes the value of one free parameter 
(eg. the initial velocity or density). Only the trans-sonic solution is stationary. 

2.4.2 Critical points of a MAES 

In a magnetized medium, there are several waves able to transport information, related 
to the different restoring forces. In a MAES, the Lorentz force couples with the plasma 
pressure gradient leading to three different MHD waves: 

• the Alfven wave (A), causing only magnetic disturbances along the unperturbed 
magnetic field Bo and of phase velocity 




• two magnetosonic waves, the slow (SM) and the fast (FM), involving both magnetic 
disturbances and plasma compression (rarefaction), of phase velocity 

vL,FM = ^ (^i + c's T ^{vl + clY-Avlcjcos^e 

where 9 is the angle between the unperturbed field Bg and the direction of propa- 
gation of the wave (the disturbance). 

In ideal MHD regime (in the jets), these three waves can freely propagate. This provides 
three singularities along each magnetic surface, whenever the plasma velocity equals one 
of these phase speeds. Thus, three regularity conditions must be specified per magnetic 
surface. 

Inside the turbulent accretion disc it is another story. There, the high level of turbu- 
lence maintains large magnetic diffusivities and viscosity. As a result, MHD waves are 
strongly damped and the number of singularities really present is not so clear. Within the 
alpha prescription of accretion discs (Shakura & Sunyaev 1973), viscosity appears only 
in the azimuthal equation of motion (angular momentum conservation). The presence 
of viscosity there "damps" the acoustic waves and there is no singularity related to this 
motion (although rotation is supersonic) . Conversely, there is no viscosity in the poloidal 
(radial and vertical) equations of motion: a singularity can therefore appear there, related 
to pure acoustic waves. As a consequence, in principle, it may be necessary to impose a 
regularity condition with respect to the accretion flow itself, if it becomes supersonic. 

3 Magnetized jets 

3.1 Commonly used equations 

As said previously, jets are in ideal MHD regime, namely — Ujji — i^*^ — 0. In this 
regime, mass and flux conservations combined with Ohm's law (30) provide 

Up = ^4 (40) 
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where rj{a) = -^fjl^PA is a constant along a magnetic surface^ and pA is the density at the 
Alfven point, where the poloidal velocity Up reaches the poloidal Alfven velocity Vap- The 
induction equation (31) becomes 

J^*(a) = J^-r?-^ , (41) 

where 17* (a) is the rotation rate of a magnetic surface (imposed by the disc, thus very 
close to the Keplerian value). Note that despite the frozen in situation, jet plasma flows 
along a magnetic surface with a total velocity u = {ri/iioP)B + 0*rej which is not 
parallel to the total magnetic field B. This is possible because field lines rotate faster 
than the ejected plasma. If the disc is rotating counter-clockwise, the field lines will be 
trailing spirals (fi* > il, ie. i?^ < with Bz > 0), while ejected plasma will rotate in 
the same direction as in the disc. Magnetic field lines and plasma stremlines are thus two 
helices of different twist. Jet angular momentum conservation simply writes 

^,r\ = nr^-^ (42) 

where va is the Alfven radius. Above the disc, the turbulent torque vanishes and only 
remains a magnetic accelerating torque. The first term on rhs is the specific angular 

momentum carried by the ejected plasma whereas the last tern can be understood as the 
angular momentum stored in the magnetic field. The total specific angular momentum 
L(a) = U^r\ is an MHD invariant. The Alfven radius can be interpretated as a 
magnetic lever arm, braking down the disc. The larger the ratio r^/ro, the larger the 
magnetic torque acting on the disc at the radius Tq. 

Hereafter, we focus only on adiabatic jets (for which thermal effects can still be non 
negligible). Usually, instead of using the other two components of the momentum conser- 
vation equation, one uses the Bernoulli equation (obtained by projecting Eq. (29) along 
the poloidal direction, ic. Bp) and the transverse field or Grad-Shafranov equation (ob- 
tained, after quite a lot of algebra, by projecting Eq. (29) in the direction perpendicular 
to a magnetic surface, ie. Va). For a jet of adiabatic index 7, Bernoulli equation writes 

E{a) = '^ + H + ^G-^*r^^ (43) 

where n is the total plasma velocity and H = (7/7 — l)P/p is the gas enthalpy. This 
equation describes the acceleration of matter along a poloidal magnetic surface, namely the 
conversion of magnetic energy and enthalpy into ordered kinetic energy. Grad-Shafranov 
equation of an adiabatic jet is 

V.(m^-l)^ = J^-J^^ + iSlr^-^^rl/^ - 
p^of I da da da 

Bl + w?Bld\ur] 



+ 



Po da 




where = u^/Vap is the Alfvenic Mach number and = jksT / /imp is the jet sound 
speed. This awful equation provides the transverse equilibrium (ie. the degree of colli- 
mation) of a magnetic surface. A simpler-to-use and equivalent version of this equation 
is 

(1 - m^)^ - P + — - pVx$G + (P^^r - -^)V^r = (45) 



^Any quantity Q verifying ■ VQ = is a constant along a poloidal magnetic field line, hence a MHD 
invariant on the corresponding magnetic surface. 
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where V_l = Va • V/| Va| provides the gradient of a quantity perpendicular to a magnetic 
surface {V±Q < for a quantity Q decreasing with increasing magnetic flux) and TZ, 
defined by 

1 _ Va {Bp • V)Bp 



n |Va| 52 



(46) 



is the local curvature radius of a particular magnetic surface. When TZ > 0, the surface is 
bent outwardly while for 7?. < 0, it bends inwardly. The first term in Eq.(45) describes the 
reaction to the other forces of both magnetic tension due to the magnetic surface (with the 
sign of the curvature radius) and inertia of matter flowing along it (hence with opposite 
sign). The other forces are the total pressure gradient, gravity (which acts to close the 
surfaces and deccelerate the flow, but whose effect is already negligible at the Alfven 
surface), and the centrifugal outward effect competing with the inwards hoop-stress due 
to the toroidal field. 

To summarize, an astrophysical jet is a bunch of axisymmetric magnetic surfaces 
a(r, z) = Const., nested one around the other at different anchoring radii Vo- The magnetic 
flux distribution a{ro) is unknown and is therefore prescribed, whereas the shape r{z) of 
the magnetic surface is self-consistently calculated. Each magnetic surface is characterized 
by 5 MHD invariants: 

- r]{a), ratio of ejected mass flux to magnetic flux; 

- ri* (a) , the rotation rate of the magnetic surface; 

- L{a) = fl*r\, the total specific angular momentum transported; 

- E(a), the total specific energy carried away; 

- K{a), related to the specific entropy S{a). 

A magnetized jet is then described by 8 unknown variables: density p, velocity u, 
magnetic field B (flux function a and toroidal fleld B^), pressure P and temperature 
T. There are 8 equations allowing us to solve the complete problem: (25), (32), (34), 
(40), (41), (42), (43) and (44) or its more physically meaningful version (45). Since 
the ejected plasma must become super-FM, 3 regularity conditions have to be imposed, 
leaving the problem with 5 free boundary conditions (at each radius To). Studies of 
disc-driven jets usually assume that matter rotates at the Keplerian speed, namely f^o = 
^k{i'o) = \/ GM/r'^ and that jets are "cold" (negligible enthalpy, K{a) = 0) or choose an 
arbitrary distribution K(a). In both cases, it leaves the problem with only 3 free and 
independent boundary conditions that must be specifled at each radius^. 

Magnetized jets are such complicated objects that only gross properties are known. For 
example, we know that a non- vanishing asymptotic current will produce a self-confinement 
of some field lines (Heyvaerts & Norman 1989). But the exact proportion of collimated 
field lines depends on "details" (transverse gradients of inner properties, outer pressure). 
The distance at which jets become collimated, the jet radius and opening angle, the 
velocity and density transverse distributions still remain to be found in full generality. 
This is the reason why there are so many different works on jet dynamics and why each 
authors use their own "relevant" parameters. Following Blandford &; Payne (1982), we 
introduce the following jet parameters 

O 7-2 Bt 



^In numerical MHD computations, one usually prescribes the density p"'"(ro), vertical velocity u+(ro) and 
magnetic field B+(ro) distributions (eg. Ouyed & Pudritz 1997, Krasnopolsky et al. 1999). In those self-similar 
jet studies where only the Alfven point has been crossed, 4 (if K{a) = 0, Blandford & Payne 1982) of 5 (if 
K{a) ^ 0, Contopoulous & Lovelace 1994) variables remain free and independent (p"*", u+, and P+). 

In both cases, the superscript stands here for the belief that this boundary condition corresponds indeed 
to the disc surface. 
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K = r}^- ~ ^2 P^z (47) 

refers to quantities evaluated at the dise midplane, "+" at the disc surface 
and "A" at the Alfven point. The first parameter, A, is a measure of the magnetic lever 
arm that brakes down the disc while n measures the ejected mass flux. Another parameter 
is usually introduced, related somehow to the jet asymptotic behaviour. We use 

^*fA ■>i/2 / .o\ 

^A = Yi ^nyl^—— (48) 

which measures the ratio of the rotational velocity to the poloidal jet velocity at the Alfven 
point. Such a parameter characterises the magnetic rotator: cold jets require loa > 1 to 
become trans- Alfvenic (Pelletier & Pudritz 1992, Ferrcira 1997, Lcry et al. 1999). Note 
also that its value depends on pure geometrical effects, namely the way the magnetic 
surface opens. As a consequence, spherical expansion of field lines is probably a special 
case, leading to particular relations between jet asymptotic behaviour and its source. 

3.2 Some aspects of cold jets physics 
3.2.1 Energetic requirements 

It is quite reasonable to assume that the magnetic energy density in the disc is much 
smaller than the gravitational energy density. As a consequence, the rotation of the disc 
drags the magnetic field lines which become then twisted. This conversion of mechanical 
to magnetic energy in the disc gives rise to an outward poloidal MHD Foynting flux 

Smhd,p = = -n,rB^^ (49) 

Mo Mo 

which feeds the jets and appears as the magnetic term in the Bernoulli integral (43). 
Another source of energy for the jet could be the enthalpy H (built up inside the disc and 
advected by the ejected flow) or another local source of heating Q (eg. coronal heating). 
We are mainly interested here in "cold" jets, where those two terms are negligible (see 
Sect. 5.3). Bernoulli equation can be rewritten as 

E{a) = f + $e// + H (50) 

where the effective potential is ^^ff = — ^-^^*r'^ with $7 = $7^,(1 — 5) and ~ Qo- 
The function g is much smaller than unity at the disc surface then increases along the 
jet (if the jet widens a lot, g ^ 1). Starting from a point located at the disc surface 
{ro,z = 0), matter follows along a magnetic surface and must move to another point 
(ro + Sro,z). This can be done only if a positive poloidal velocity is indeed developed. 

Making a Taylor expansion of $e//) one gets ^ Hg — H + ^{Sdr^ — z^) > 0, which 
translates into the condition 

ta.. = ^<v/3(l + ?^)"\ (51) 

Thus, cold jets (negligible enthalpy) require fleld lines bent by more than 30" with respect 
to the vertical axis at the disc surface (Blandford & Fayne 1982). The presence of a 
significant enthalpy {Hg large) is obviously required if this condition is not satisfied. 
Bernoulli equation (43) gives a total energy feeding a cold jet 

E{a) = ny, + n^iny^ - n^rl) = ^(2A - 3) (52) 
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which is directly controlled by the magnetic lever arm A. Cold, super-FM jets require 
therefore A > 3/2. If all available energy is converted into kinetic energy, ejected plasma 
reaches an asymptotic poloidal velocity Urx>{a) = \/2E{a). 

3.2.2 Relevant forces and current distributions 

So, rotation of open field lines produces a shear {B'^ < 0) that results in an outward flux 
of energy. Once matter is loaded onto these field lines, it will be flung out whenever there 
are forces overcoming the gravitational attraction. Since matter flows along magnetic 
surfaces, one must look at the projection of all forces along these surfaces. We obtain for 
the Lorentz force, 

Bn 

F|| = -;^V||/ (53) 
II 27rr " 

where V|| = V • Bp/ Bp and / = 2iTrB^/iJLo is the total poloidal current flowing inside 
the magnetic surface. Hence, jets are magnetically-driven whenever V||/ > is fulfilled, 
namely when current is leaking through this surface. This quite obscur condition de- 
scribes the fact that magnetic energy is being converted into kinetic energy: the field lines 
accelerate matter both azimuthaly {F^ > 0) and along the magnetic surface (Fy > 0). 
The difference with the Barlow wheel lies in the possibility to convert magnetic energy 
into jet (bulk) kinetic energy. 

The jet transverse equilibrium depends on the subtle interplay between several forces 
(see Eq. (45)). The transverse projection of the magnetic force, namely 

F^ = BpJ^-^V^I (54) 

where V_l = (Va • V)/|Va|, shows that it depends on the transverse current distribution. 
Thus, the degree of jet collimation (as well as plasma acceleration) depends on the overall 
electric current circuit. Any bias introduced on this circuit can produce an artificial force 
and modify diagnostics on jet collimation. 

Eventually matter becomes no longer magnetically accelerated and Vy / = is satisfied. 
This implies two possible asymptotic current distributions. If jets are force-free (ie. Jp and 
Bp parallel), then there is a non- vanishing asymptotic current providing a self-collimating 
pinch (and one must worry about how the electric circuit is closed). Or, jets become 
asymptotically current-free {I^ = 0) and another cause must then be responsible for 
their collimation (either inertial or external pressure confinement). This last alternative 
has something appealing for magnetic fields would then have a major influence only at 
the jet basis, becoming dynamically negligible asymptotically (cf Sect. 1.1). 

3.3 Numerical simulations: what can be learned of them ? 

There have been a lot of numerical studies of MHD jet propagation and their associated 
instabilities. Here, I focus only on those addressing the problem of jet formation from 
accretion discs. Although some attempts have been made to model accretion discs driving 
jets, difficulties are such that nothing realistic has been obtained yet (Shibata k. Uchida 
1985, Stone & Norman 1994, Kudoh et al. 1998). Ejection is indeed observed, but no one 
can tell whether these events are just transients or if they indeed represent some realistic 
situation. 

Another philosophy is to treat the disc as a boundary condition and, starting from an 
(almost) arbitrary initial condition, wait until the system converges towards a stationary 
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flow (Ustyugova et al. 1999, 2000, Ouyed & Pudritz 1997, 1999, Krasnopolsky et al. 
1999). Now, it is of no wonder that jets are indeed obtained with these simulations. 
Matter, forced to flow along open magnetic field lines, is continuously injected (at a rate 
p'^uf) at the bottom of the computational box. As a response, the field lines twist 
(ie. B(j) increases) until there is enough magnetic energy to propell it. If the code is 
robust enough, a steady-state situation is eventually reached, actually reproducing most 
of the results obtained with self-similar calculations (Krasnopolsky et al. 1999). Ouyed 
& Pudritz (1997) found a parameter region where unsteady solutions are produced, with 
a "knot generator" whose location seems to remain fixed. Their results may be related to 
the characteristic rccollimation configuration featured by self-similar solutions (Blandford 
& Payne 1982, Contopoulos & Lovelace 1994, Ferrcira 1997). 

Nevertheless, since the time step for computation depends on the faster waves (Courant 
condition) whose phase speed varies like it appears that no current MHD code can 

cope with tiny mass fluxes. Thus, code convergence itself introduces a bias in the mass 
flux of numerical jets. These are always very "heavy", with a magnetic lever arm r^/r^ 
not even reaching 4 (Ouyed &: Pudritz 1999, Krasnopolsky et al. 1999, Ustyugova et 
al. 1999) while cold self-similar studies obtained much larger values (with a minimum 
value around 10). This is not a limitation imposed by physics but of computers and will 
certainly be overcome in the future. 

On the other hand, the question of jet collimation is clearly one that can be addressed 
by simulations. However, the boundary conditions imposed at the box, as well as the 
shape of the computational domain, can introduce artificial forces leading to unsteady 
jets or spurious collimation (see the nice paper of Ustyugova et al. 1999). 

Anyway, the following crucial question remains to be addressed: how (and how 
much) matter is loaded into the field lines ? Or another way to put it: how is 
matter steadily deflected from its radial motion (accretion) to a vertical one (ejection) ? 
To answer this question one must treat the accretion disc in a self-consistent way. 



4 Magnetized accretion discs driving jets 

4.1 Physical processes in quasi-Keplerian discs 
4.1.1 Turbulent, Keplerian discs 

As said previously, we focus our study on Keplerian accretion discs. Such a restriction 
imposes negligible radial plasma pressure gradient and magnetic tension. We define the 
local vertical scale height as Pg = Po^^h^ where Pg is the total plasma pressure. Looking 
at the disc radial equilibrium, a Keplerian rotation rate is indeed obtained whenever the 
disc aspect ratio 

e = ^ (55) 
r 

is smaller than unity. Hereafter, we use e < 1 as a free parameter and we will check a 
posteriori the thin disc approximation (Sect. 5.4.1). 

Steady accretion requires a turbulent magnetic diffusivity for matter must cross field 
lines while accreting and rotating. Since rotation is much faster than accretion, there 
must be a higher dissipation of toroidal field than poloidal one. A priori, this implies 
a possible anisotropy of the magnetic diffusivities associated with these two directions, 
poloidal and toroidal f^. Besides, such a turbulence might also provide a radial 
transport of angular momentum, hence an anomalous viscosity To summary, at least 
three anomalous transport coefficients are necessary to describe a stationary MAES. We 
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will use the following dimensionless parameters defined at the disc midplane: 

level of turbulence 

degree of anisotropy (56) 
magnetic Prandtl number 

where va = Bo/ ^Jp^oPo is the Alfven speed. Our conventional view of 3D turbulence 
would translate into < 1, Xm ~ 1 and Vm ~ 1- But as stated before, the amount 
of current dissipation may be much higher in the toroidal direction, leading to <C 1. 
Moreover, it is not obvious that a^i must necessarily be much smaller than unity. Indeed, 
stationarity requires that the time scale for a magnetic perturbation to propagate in the 
vertical direction, h/vA, is longer than the dissipation time scale, h? /um. This roughly 
translates into > 1- Thus, we must be cautious and will freely scan the parameter 
space defined by these turbulence parameters. 
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4.1.2 How is accretion achieved ? 



The disc being turbulent, accretion of matter bends the poloidal field lines whose steady 
configuration is provided by Ohm's law (30) . At the disc midplane, this equation provides 



TZ = 



ruo 



(57) 



z=0 



where 71^ is the magnetic Reynolds number related to the radial motion Uo and l{r) is 
the characteristic scale height of the magnetic flux variation. Once the disc turbulence 
properties are given, we just need to know the accretion velocity Uo- This is provided by 
the angular momentum equation, namely 



r 



+ (V • T) • 



(58) 



where is the magnetic torque due to the large scale magnetic field ("jet" torque) and 
(V • T) • is the "viscous" -like torque of turbulent origin (possibly due to the presence of 
a small scale magnetic field) . Such an equation can be put into the following conservative 
form 



V 



rB^ 
Po 



-Br, 







(59) 



where r (V • T) • = V • rTy . Although modelling the jet torque is quite straightforward, it 
is not the case of the turbulent torque Ty and we use a prescription analogous to Shakura 
& Sunyaev (1973). Defining 



A 



jet torque 



turbulent torque 



the disc angular momentum conservation becomes at the disc midplane 



1 + A ~ 



TZr, 
Vr, 



(60) 



(61) 



Taking the conventional value Vm ~ 1, one sees that a "standard" accretion disc, which 
is dominated by the viscous torque (A ^ 1) requires straigth poloidal field lines (TZm ~ 
1, Heyvaerts et al. 1996). On the contrary, cold jets carrying away all disc angular 
momentum (A S> 1) are produced with bent magnetic surfaces so that TZm ~ A ~ 
(Ferreira & Pelletier 1995). 
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Figure 1: Sketch of a Magnetized Accretion Ejection Structure (Courtesy of Fabien Casse). 

4.1.3 How is matter defieted from accretion ? 

The poloidal components of the momentum conservation equation write 



At the disc midplane, a total (magnetic + viscous) negative torque provides an angular 
velocity slightly smaller than the Keplerian one = y/GM/r^, thereby producing an 
accretion motion (ur < 0). But at the disc surface, one gets an outwardly directed flow 
{u^ > 0) because both the magnetic tension Fr (whose effect is enhanced by the fall in 
density) and the centrifugal force overcome gravitation (J7+ > i^x)- Note that jets are 
magnetically-driven, the centrifugal force resulting directly from the positive magnetic 
torque (F+ > 0). This can be understood with a simple geometrical argument: matter 
has been loaded onto field lines that are anchored at inner radii and are thus rotating 
faster. 

But again, we assumed that matter is being loaded from the underlying layers with 
> 0. The physical mechanism is hidden in Eq. (63): the only force that can always 
counteract both gravity and magnetic compression is the plasma pressure gradient (Fer- 
reira & Pelletier 1995). Inside the disc, a quasi-MHS equilibrium is achieved, with matter 
slowly falling down (uz < 0) while accreting. However, plasma coming from an outer 
disc region eventually reaches the upper layers at inner radii. There, the plasma pressure 
gradient slightly wins and gently lifts matter up {uf > 0), at an altitude which depends 
on the local disc energetics. 

4.1.4 How is steady ejection obtained ? 

Is ejection unavoidable once all above ingredients are met ? The answer is "yes", but 
steady ejection requires another condition. 

''Namely, rotation, open field lines and some amount of diffusion allowing loading of matter. 




(62) 



(63) 
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While accretion is characterized by a negative azimuthal component of the Lorentz 
force magnetic acceleration occurring in jets requires a positive F^. Since = 
JzBf — Jf-Bz, the transition between these two situations depends mainly on the vertical 
profile of the radial current Jr = —iig^dB^/ dz, that is, on the rate of change of the 
magnetic shear with altitude. In order to switch from accretion to ejection, must 
vertically decrease on a disc scale height. This crucial issue is controlled by Eq. (31), that 
provides 

rl^Jr-rljo + r fdzEp-Vn - B^u, (64) 
Jo 

The first term on the rhs describes the current due to the electromotive force, the second 
is the effect of the disc differential rotation and the third is the advection effect, only 
relevant at the disc surface layers. Thus, the vertical profile of Jr is mainly controlled by 
the ratio of the differential rotation effect over the induced current (Ferreira & Pelletier 
1995). No jet would be produced without differential rotation for it is the only cause of 
the vertical decrease of Jr- However, the counter current due to the differential rotation 
cannot be much bigger than the induced current in the disc, otherwise Jr would become 
strongly negative and lead to an unphysically positive toroidal field at the disc surface. 
Thus, steady state ejection is achieved only when these two effects are comparable, which 
translates into 

For A > Ac, matter is spun down at the disc surface, while for A < Ac there is not 
enough energy to propell the large amount of mass trying to escape from the disc. Thus, 
equation (65) is a necessary condition for stationarity. 



4.1.5 From resistive discs to ideal MHD jets 

As matter is expelled off the disc (by the plasma pressure gradient) with an angle Oy,^ = 
arctan(n,./uz), magnetic stresses make it gradually flow along a magnetic surface (with 
Up II Bp). Indeed, Ohm's law (30) can be written 

(tan^R \ 
" V ^^^^ 

where Ob^ = &Tcian{Br/ B^)- As long as < Osp, the toroidal current remains positive. 
This maintains a negative vertical Lorentz force (that decreases u^) and a positive radial 
Lorentz force (that, along with the centrifugal term, increases Ur), thus increasing 9up- 
If 9up > 9 Bp-, the toroidal current becomes negative, lowering both components of the 
poloidal Lorentz force and, hence, decreasing 9up- Therefore, in addition to the vertical 
decrease of the magnetic diffusivity (for its origin lies in a turbulence triggered inside the 
disc), there is a natural mechanism that allows a smooth transition between resistive to 
ideal MHD regimes. 



4.2 Disc-jets interrelations 
4.2.1 Dimensionless parameters 

The accretion disc is defined by 11 variables, ie. the same 8 as in the ideal MHD jets plus 
the 3 transport coefficients. All these quantities can be calculated from their values at 
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the disc equatorial plane. At a particular radius they write 



rp _ GMmp_ 2 



dz 



:=0 
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n — Mao 



Bo 

dz 
J, 



1/2 



,0 



where = 2g/x^^ = aye{l + A) 



where 5 



where 
where 
where 



^ + 2(2 - /3) + nTZr, 



1/2 



q = 



2/.1/2 ^ 
= Vm{l + A) 



where i^^ = a^^Kh'^-, = V^M/r'^, Mao = ^a{ro), /3 = dlna/d\nro provides the mag- 
netic flux distribution and A is constrained by Eq. (65). Thus, there are 3 more parameters 
in addition to the previous 4 ones (e, am, Vm, Xm)'- the magnetic flux distribution /3, the 
magnetic held strength and the ejection index ^. The 3 regularity conditions arising at 
the 3 critical points met by the ejected plasma along each magnetic surface provide the 
value of 3 disc parameters (precisely, their values at the anchoring radius Tq). Inside our 
cold approximation, we therefore expect to fix the values of jx and ^ as functions^ of 
the free parameters e, a^,, Vm and Xm- As a consequence, jet properties (ie. invariants 
as well as the asymptotic behaviour) arise as by-products of these parameters. 

Using the set of ideal MHD equations, mass conservation gives a mass flux leaving 
the disc surface p'^u^ — £,epoUo related by ^ to the accretion mass flux. Then, angular 
momentum conservation provides the following exact relations 
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Both magnetic lever arm A and mass load n are therefore determined by the ejection index 
^ (which is directly related to the exact value of the toroidal fleld B'^ at the disc surface^). 
One way to understand this point is to compute the ratio a of the MHD Poynting flux to 
the kinetic energy flux. At the disc surfaces this ratio writes 
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The ejection index appears to be also a measure of the power feeding the jets. Unless 
^ is of order unity, the magnetic field completely dominates matter at the disc surface 
(Eq. (65) forbids 1 >> A ~ ^). 



*An important remark. The field strength cannot be too large (/i >> 1), otherwise there will be no vertical 
equilibrium possible. On the other hand, if it is too small (/k << 1), then the disc is prone to the magneto- 
rotational instability (Balbus & Hawley 1991). Therefore, we expect ~ 1 in steady-state MAES. However, 
the allowed values of ^ strongly depend on the (subtle) vertical equilibrium and its interplay with the induction 
equation. 

®The estimate B'^ ~ —qBo is roughly acurate (by less than a factor 2) but the jet asymptotic structure 
highly depends on its precise value (Ferreira 1997). 
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4.2.2 Constraints on the ejection index of cold MAES 



Can we provide any general constraint on the allowed range of ^ ? This is the trickiest 
question about MAES and here follows some analytical arguments in the simple case of 
cold jets from isothermal discs (Ferreira 1997). 

The maximum value £,max is constrained by the jet capability to accelerate ejected 
matter, hence momentum conservation. Let's consider two extreme cases. Imagine there 
is such a tiny fraction of mass ejected that it takes almost no energy to accelerate it 
until the Alfvcn point. But as some acceleration has been provided anyway, we can write 

< <7+. On the other extreme, a huge amount of matter is expelled off the disc, which 
has hardly enough energy. In this most extreme case, no more energy is left after the 
Alfven point and m^^j, = 1. Gathering these two conditions provides the constraint 

Tr^<-A<Y^ (69) 

which shows that cold jets (1) require a;^ > 1 (fast rotators) and (2) display a maximum 
ejection efHciency of ^max = — 0.15. Higher ejection indices are inconsistent with 

steady-state, trans- Alfvcnic, cold jets. 

The minimum value ^min of the ejection index is constrained by the disc vertical 
equilibrium. The diminishing of the ejection efficiency ^ is obtained by increasing the 
magnetic compression, especially through the radial component (via A, increasing as ^ 
decreases). Now, to maintain the vertical balance while bending the field lines, but without 
increasing the plasma pressure, one must decrease the magnetic field amplitude (parameter 
jj). But then, there is a non- linear feedback on the toroidal field induction. Indeed, 
as jj, decreases, the effect due to the differential rotation decreases also, leading to an 
increase of the toroidal field at the disc surface. This causes an increase of the toroidal 
magnetic pressure, hence an even greater magnetic squeezing of the disc. Below ^mmj 
no vertical equilibrium is possible. Providing a quantitative analytical expression of how 
much matter can actually be ejected, ie. the value of ^rnm; is out of range. Indeed, in 
the resistive upper disc layers, all dynamical terms are comparable in Eq. (63). A 
careful treatment of the disc vertical balance is therefore badly needed. Thus, finding 
the correct parameter space of a MAES forbids the use of crude approximations, like 
puz = constant (Wardle &; Konigl 1993), using strict hydrostatic balance (Li 1995) or any 
other prescription mimicking the induction equation (Li 1996). 



4.2.3 Jet asymptotic structure 

Can we relate the jet asymptotic structure to the disc parameters without solving the full 
set of MHD equations ? Naively, one would say that the larger A the larger jet asymptotic 
radius, if some cylindrical collimation is achieved. Next section, we will see that this last 
issue is far from being obvioTis. Anyway, we can still safely say that if some current is still 
available after the Alfven point, then magnetic acceleration will probably occur, along 
with an opening of the magnetic surfaces. In fact, the larger B^, the larger B^^^. This, 
in turn, ensures that jets will provide a big acceleration and open up a lot. For a cold jet, 
the ratio Ia/I+ of the remaining current I a to the current provided at the disc surface 
/+ is 

Ia vaB^^a , 2 -, 3 1 2 

T~ = ^ - 9A where ^^ = 1- - + 70 

1+ ToBJ A uj% A-^/^(l + 4/r^)V2 

The expression of qa is the Bernoulli equation evaluated at the Alfven point and 

^^^^^3/2^m(^d^ 
sm cpA 
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where (j)A is the local angle between the Alfven surface (defined hy z = 2;a(^a)) and the 
vertical axis, and 6a the opening angle estimated at the Alfven point. Both angles are 
determined by the resolution of the Grad-Shafranov equation, which takes into account 
radial boundary conditions. Note that this expression is only valid for a conical Alfven 
surface, a geometry which arises naturally when jet parameters vary slowly across the jet. 
Since cold jets require fast magnetic rotators (loa > !)> a necessary condition for trans- 
Alfvenic jets is kX^/'^ > 1, which translates into A > (Aai,e)~^. Thus, jets launched from 
discs with a dominant viscous torque A < 1 require huge magnetic lever arms, namely 
A S> £~^. This is most probably forbidden by the disc vertical equilibrium that would 
not survive such a strong magnetic pinching. So, cold disc-driven jets are presumably 
carrying away a significant fraction of the disc angular momentum (A > 1). 

4.2.4 Summary 

From the preceeding general analysis, we can a priori expect two extreme cold configu- 
rations from quasi-Keplerian discs: 

Type I, where large toroidal currents at the disc midplane correspond to large mag- 
netic Reynolds numbers TZm ~ » 1 and a dominant magnetic torque A ~ £~^. 
This configuration would be achieved for isotropic turbulence, Vm ~ 1 and Xm ^ 1- 

Type II, where the dominant source of toroidal currents is at the disc surfaces, corre- 
sponding to straight field lines inside the disc {TZm ~ 1) that become bent only 
at the surface. These surface currents come from the electromotive force {J^ ~ 
— {urBz/rjm)'^), due to the presence of a large viscous torque (A ~ 1) allowing a non 
zero accretion velocity at the disc surface. Such a configuration would be achieved 
for an anisotropic turbulence, Vm > 1 and Xm ~ £• 

At this stage, I hope the reader has achieved an understanding of the relevant physical 
mechanisms inside a keplerian accretion disc driving jets (question 1, Sect. 1.4). The disc 
physical conditions (question 2) are described by the MAES parameter space and thus, 
require the treatment of the complete set of MHD equations. As a "by-product", we will 
hopefully have the answer of the last question (jet properties). 

5 Self-similar solutions of MAES 

5.1 Self-similar Ansatz and numerical procedure 

Solving the full set of MHD equations requires heavy 2D or 3D numerical simulations. 
However, looking for special solutions will allow us to transform the set of partial dif- 
ferential equations (PDE) into two sets of ordinary differential equations (ODE) with 
singularities. Gravity is expected to be the leading energy source and force in accretion 
discs. Thus, if MAES are settled on a wide range of disc radii, magnetic energy den- 
sity probably follows the radial scaling imposed by the gravitational energy density. The 
gravitational potential writes in cylindrical coordinates 



Since the disc is a system subjected to the dominant action of gravity, any physical 
quantity A{r, z) will follow the same scaling, namely A{r, z) = GA{r)fA{^)- Since gravity 
is a power law of the disc radius, we use the following self-similar Ansatz 




(72) 




(73) 
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where x = z/h{r) = z/er is our self-similar variable and is the MAES outer radius. 
Because all quantities have power law dependencies, the resolution of the "radial" set 
of equations is trivial and provides algebraic relations between all exponents. The most 
general set of radial exponents allowing to take into account all terms in the dynamical 
equations (ie. no energy equation) is: 

^ = i + | "'=(-1 

otBr = OiB^ = OiB^ = f3 - 2 ap = ap - 1 

_ _ _ _1 _ _ _ 1 

As an illustration, the solutions obtained by Blandford &: Payne used (3 = 3/4, ie ^ = 0. 
Note also that the disc scale height must verify h{r) = er. Such a behaviour stems only 
from dynamical considerations, ie. the vertical equilibrium between gravity and mag- 
netic compressions and plasma pressure gradient. However, the energy equation provides 
another constraint that is usually incompatible with such a scaling (see eg. Ferreira & 
Pelletier 1993). We will come back to this issue later on. 

All quantities fA{x) are obtained by solving a system of ODE which can be put into 
the form 
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V 



/ 1 \ 



\ dx / 



where M is a 8x8 matrix in resistive MHD regime, 6x6 in ideal MHD (see Ferreira &: Pel- 
letier 1995). A solution is therefore available whenever the matrix M is inversible, namely 
its determinant is non-zero. Starting in resistive MHD regime, det M = whenever 

V'^iV'^ -Cl) = Q , (74) 

where Cg is the sound speed and V = u-n\s the critical velocity. The vector 

_ &z ~ xecr , . 

n = -77-. — TT^ (75) 

provides the direction of propagation of waves that are consistent with our axisymmetric, 
self-similar description. Therefore, close to the disc, the critical velocity is y ~ n^, whereas 
far from the disc it becomes V ~ Uj. (no critical point in the azimuthal direction). Thus, 
inside the resistive disc, the anonalous magnetic resistivities produce such a dissipation 
(presence of high order derivatives) that the magnetic force does not act as a restoring 
force and the only relevant waves are sonic. Note also that the equatorial plane where 
y = is also a critical point (of nodal type since all the solutions must pass through it). 
This introduces a small difficulty, since one must then begin the integration slightly above 
it. In the ideal MHD region, det M = whenever 

iV' - V^mW^ - V^mW - VLf = (76) 

namely, where the flow velocity V successively reaches the three phase speeds Vsm, Vati 
and Vfm, corresponding respectively to the slow magnetosonic wave, the Alfven wave and 
the fast magnetosonic wave. The phase speeds of the two magnetosonic modes have the 
usual expression, namely = \ (c"^ + Vj^ =f \J (C| + _ AC^Vj^) where V^t 

is the total Alfven speed and Vau = ^Ap ■ n. Note that the condition V = Vau is equivalent 
to Up = Vap, which shows that this is the usual Alfvenic critical point encountered in jet 
theories. It is also noteworthy to remark that the multiplicity of this root implies that at 
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the Alfvenic point, both first and second order derivatives of the physical quantities are 
imposed by the regularity condition. 

How do we proceed ? Starting slightly above the disc midplane where all quantities are 
known, we propagate the resistive set of equations using a Runge-Kutta (better a Stoer- 
Burlisch) solver. We do this for fixed values of the four free parameters {£,am,Xm,'Pm) 
and some guesses for fj, and ^. As x increases, the flow reaches an ideal MHD regime 
and we shift to the corresponding set of equations. Care must be taken in order to not 
introduce jumps in the solution while doing this. If, for the chosen ejection efficiency ^, the 
value fj, of the magnetic field is too large, the overwhelming magnetic squeezing leads to a 
decrease of u^- If, on the contrary, /x is too small, the plasma pressure gradient becomes 
far too efficient and leads to infinite vertical acceleration. Using these two criteria and 
fine-tunning we can approach the SM point so close that a simple linear extrapolation 
on all quantities allow us to safely cross the singularity. This leapfrog must not introduce 
any discontinuity (to some tolerance) on all jet invariants. By doing so, we obtain a trans- 
SM solution. This is done for a chosen value of ^ that may not be the critical allowing a 
trans- A solution. If ^ < the magnetic tension wins over the outwardly directed tension 
produced by the ejected, rotating plasma. As a result, the magnetic surface closes which 
leads to a decceleration and the Alfven point is not reached. On the other hand, ^ > 
produces an over-opening of the magnetic surface leading to the unphysical situation where 

goes to zero. Again, once close to the Alfven point (typically, m = 1 by 1%), we do 
a leapfrog. Thus, fine-tunning ^ and, for each guess, finding the critical value of n, allow 
us to obtain trans-SM and super-Alfvenic solutions. No trans-FM solution connected to 
the disc has been found yet. 

5.2 "Cold" configurations 

A cold configuration is defined by a negligible enthalpy at the disc surface (at each mag- 
netic surface). Since discs are quasi-keplerian, cold jets are obtained with an isothermal 
(Wardle k Konigl 1993, Ferreira & Pelletier 1995, Li 1995, 1996, Ferreira 1997) or adia- 
batic (Casse & Ferreira 2000a) vertical profile frix) of the temperature. 

5.2.1 General behaviour 

The general behaviour is exactly that expected in the disc. Both magnetic and (if non 
negligible) turbulent viscous torques extract energy and angular momentum from disc 
plasma. The accretion velocity depends on the amount of magnetic diffusivity. As plasma 
is being accreted from the outer disc regions, it is slightly converging towards the disc 
midplane (because of both tidal and magnetic compressions) . Matter located (locally) at 
the disc surface feels a strong outwardly directed magnetic tension and a positive azimuthal 
torque, both arising from current consommation (V||/ > 0): accretion is stopped and 
reversed. More or less simultaneously, a positive vertical velocity is provided by the 
plasma pressure gradient. 

Ejected matter leaves the disc with a vertical velocity initially much smaller than the 
local sound speed, ~ mg^Cs- It gets however very quickly accelerated as it leaves the 
resistive MHD zone. The SM point lies typically between 1 and 2 scale heights, usually 
at the very beginning of the ideal MHD zone. Until the Alfven point, plasma is almost 
co-rotating with the magnetic field lines, behaving like a rigid funnel. The Alfven point is 
far away above the disc, at an altitude za ^ ta = X^^'^ro ^ h{ro). After its crossing, there 
is a sudden opening of the magnetic surfaces. This is due to the centrifugal force which is 
now enhanced by the tension provided by the super-Alfvenic flow {m? > 1, see Eq. 45). 
This opening of the magnetic surfaces is controlled by the quantity ua'- the larger oja, the 
larger maximum jet radius. Or, the less current used in the sub-A region and the more 
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Figure 2: Left: Components of the jet poloidal velocity Up and ratio of the poloidal to the 
toroidal velocities, measured along a magnetic surface and different ejection indices: ^ = 0.005 
(sohd hne), ^ = 0.01 (dotted line), ^ — 0.02 (short-dashed line) and ^ = 0.05 (long-dashed 
line). Note that the poloidal velocity is almost zero at the disc surface (a; ~ 1). Right: Ratio a 
of the MHD Poynting flux to the kinetic energy flux, logarithm of the Alfvenic Mach number 
m^, total poloidal current / embraced by the magnetic surface (normalized to its value at the 
disc surface) and plasma density logarithm along a magnetic surface. 



remains in the super-A region. 

As the magnetic surface opens up, plasma drags along the field lines which produces a 
large toroidal field. This is the consequence of ideal MHD since matter is rotating slower 
than the field lines (il < Q*)- As the toroidal field increases, the "hoop stress" becomes 
more and more important until it overcomes the centrifugal force. This is unavoidable if 
the jet can open up freely in space (negligible outer pressure) and if its inner transverse 
equilibrium allow it: the centrifugal effect decreases faster with radius than the hoop- 
stress. As a result, the cold jet asymptotic transverse equilibrium is simply given by 

/ r2 \ ^2 

V 2//o ) HoT 

where Pext must be understood here as the pressure of the material located on the axis or 
outside the jet. A strict cylindrical collimation requires therefore a perfect (and fragile) 
matching between the hoop-stress and the total pressure gradient. In our self-similar 
solutions, Pf,xt = Ei'iid the inner magnetic field gradient is not enough to balance the 
hoop-stress. As a consequence, all solutions recoUimate (refocus) towards the jet axis. 

This seems to be a (quite) general result of MHD jets that can freely expand in space. 
For example, all Blandford & Payne disc wind solutions and X-wind solutions obtained 
by Shu et al. (1995) also recollimate^'^ . However, Contopoulos & Lovelace (1994) as well 
as Ostriker (1997) obtained solutions within the same self-similar ansatz that did not 
recollimate: recollimation is therefore not a feature of self-similarity alone. On the other 
hand, Pelletier & Pudritz (1992) found recollimating solutions that are not self-similar. 

^°Shu et al. imposed however a cylindrical asymptotic collimation by assuming an equilibrium between the 
jet hoop-stress and an inner magnetic pressure, provided by a poloidal field located on the axis. 
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Figure 3: Left: Poloidal magnetic surfaces for e = 0.01 and = 1 and different ejection 
indices (fience cua)- $ ^ 0.05, 0.04, 0.03, 0.02, 0.012, 0.01, 0.009, 0.007 and 0.005 (the maximum 
radius increases with decreasing ^). The thick line connects the Alfven points of each surface, 
anchored at a radius Vo- Right: Critical speed V = u ■ n and phase speeds Vsm, Vati, Vfm 
(corresponding respectively to the SM, Alfven and FM waves) along a magnetic surface. All 
speeds are normalized to the disc sound speed flRhp- The SM point lies just above the disc 
surface at zsm ~ h, while the Alfven point is at za ~ ta (here ~ lOro). The solution shown 
here does not cross the last (FM) critical point (although Up > Vfm)- 



In fact, it is possible to show that if the conditions 

dlnX ^ ^ dln{pA/po) ^ dlnpo ^ ^ 

dlnro din To d In ro ^ 

are satisfied in a cold, disc-driven jet, then recollimation would have been impossible 
(Ferreira 1997). Evidently; such conditions are violated in self-similar solutions: A and 
Pa/Po remain constant throughout the jet. However, this analytical analysis shows that 
non self-similar jets produced from a large range of disc radii and displaying no strong 
gradients of these quantities may be prone to recollimation (as in Pelletier &; Pudritz). 

What happens to our solutions after they recollimate ? When the jet transverse 
equilibrium enforces recollimation, almost all available angular momentum has already 
been transfered to matter and O ~ 0*r^/r^. This implies that the current also goes to 
zero 1^0 (and not a constant) and so does the toroidal field. The recollimating matter 
drags the field lines along with it, severely reducing the pitch of the magnetic helix. 
Before Bfj) reaches zero, our self-similar solution meets the last critical point, namely 
V ^ Ur = Vfm ~ VAt (note that Up > Vfm already before recollimation). Vlahakis et al. 
(2000) obtained recently super-FM solutions by playing with the location of the Alfven 
surface and the jet polytropic index. But these solutions are terminated in the same way 
as those displayed here. Such a behaviour remains unexplained and may be due to the 
self-similar form of the solutions. Anyway, such super-FM solutions could safely produce 
an oblique shock^^ leading to a time-dependent readjustment of the whole structure. This 
will not alter the underlying steady-state solution, for no signal can propagate upstream. 



^^As suggested by Gomez de Castro & Pudritz (1993) and maybe seen in numerical simulations of Ouyed & 
Pudritz (1997). 
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Figure 4: Parameter space of cold, adiabatic MAES, for am = 1 and e = 0.1. The shaded 
areas correspond to the location where numerical solutions could be found. Top: thick dashed 
lines show theoretical limits, super-SM (upper) and super-A (lower) conditions. Levels for the 
corresponding jet parameters n and cua are also displayed. Bottom: the corresponding jet k — \ 
parameter space. 

Although Vlahakis et al. solutions are not connected with the disc, their work gives an 
indication that introducing another degree of freedom (the polytropic mdcx value) might 
indeed provide trans-FM solutions. But in any case, at this stage, a shock seems the 
unavoidable fate of this mathematical class of solutions. 

5.2.2 Parameter space of cold, adiabatic MAES 

The parameter space of cold MAES is obtained by varying the set of free disc parameters 
{e,am,'Pm,Xm)- We choose to fix the values of both e and am and represent the param- 
eter space with the remaining parameters (Fig. 4). Numerical solutions are only found 
inside the shaded areas, where we also plot levels of two jet parameters cua and n. This 
region is embedded inside a larger region (thick dashed lines), obtained by two analytical 
constraints. The first one arises from the requirement that jets become super-SM and is 
thus related to the disc vertical equilibrium. The second emerges from the requirement 
that jets must become super- Alfvcnic, namely cua > 1, providing the lower limits in the 
same plots. These two analytical constraints strongly depend on both am and s. Cold, 
adiabatic MAES have the following properties: 

• The parameter space is very narrow, with typical values ^ ~ 10~^ and /i ~ 1. No 
solution has been found outside the range 0.001 < e < 0.3 and 0.3 < a^n < 3. The 
corresponding jet parameters lie in the range 10 < A < 100 and 0.001 < k < 0.03. 

• The parameter space shrinks considerably with am, because of the extreme sensibility 
of to it (Casse & Ferreira 2000a). Note however that a„ = amlJ^^'^Vm is usually smaller 
than am- 
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• No solution has been found with a dominant viscous torque (A <C 1): the reason 
Hes in the imposed geometry of the Alfvcn surface (conical). On the other hand, high-cj^ 
solutions (those with large jet radius) exist only for magnetically-dominated discs (A ^ 1). 

• The required turbulence anisotropy increases with both A and am, following the 
scaling Xm ~ Aea^a^ provided by Eq. (65). 



5.3 "Warm" configurations 

5.3.1 Entropy generation inside the disc 

As seen in Sect. 2.3.2, the disc energy equation (33) is such a mess that all works on 
accretion discs used simpified assumptions. In fact using an isothermal or adiabatic pre- 
scription for the temperature vertical profile may be an over simplification, that may have 
led to such a small parameter space. One way to tackle the energy equation is to solve 

pT— = pTup ■VS = Q (79) 

where the entropy source Q = F — A is prescribed and describes the local net effect of all 
possible heating T and cooling A terms. The sources of heating are: 

^ Tgff = rjmJ^ + v'm'^p + "Hv \r'VO,\'^ the effective Joule and viscous dissipations; 

^ Tturb due to turbulent energy deposition, not described by anomalous coef- 
ficients; 

^ Text some external source of energy (like protostellar UV or X-rays, or cosmic 

rays) . 

On the other hand, the cooling sources are 

^ ^rad = V • Srad radiative losses in optically thick or thin media; 
^ ^turb turbulent transport that may be described by kinetic theory or due to 
large scale motions like convection (in Eq. (33) V • qturb = ^tm-b ~ '^turb)- 

In our simplified approach, we prescribe both the shape and amplitude of Q, but 
consistently with energy conservation. We make therefore two assumptions: (1) there 
is no net input of turbulent energy in the volume V occupied by the MAES, namely 
Pturb = Iv ~ ^)turb '^^^ ~ ^ (2) the power deposited by any external medium is 
negligible, ie. P^xt = /v ^extd'^'^ = 0. With these assumptions, the only remaining source 
of energy is accretion of the laminar flow: turbulence can only redistribute energy from 
one place to another one. 

Accretion is possible because of the torque due to the mean magnetic field (jet) and the 
turbulent ("viscous") torque. Accretion energy released by the first torque is converted 
into a MHD Poynting flux leaving the disc (and feeding the jets) and heat through local 
Joule dissipation. Accretion energy released by the "viscous" torque is conventionaly 
thought as being converted into heat through dissipation. But note that if such a torque 
arises from field lines connecting two disc radii, one would then expect also an outward 
flux of energy, which would be dissipated above the disc surface (Heyvaerts &; Priest 1989, 
Miller Sz Stone 2000, Machida et al. 2000). Anyway, these two heating sources (Tg//) 
build up a local gas thermal energy reservoir which decreases because of the local cooling 
terms. Thus, the total power related to dissipation inside the disc (ic. not directly put into 
the jets) is Pdiss = /y ^effd^"^- In a conventional picture of accretion discs, such a power 
is finally radiated away, either at the disc surfaces only, or also in some chromosphere. 
Here, we assume that a fraction / of this power is in fact not lost by the plasma but 
provides an extra source of entropy Q, namely 
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Figure 5: Left: Variation of jet parameters with /. As / increases, jets get hotter, denser 
and with a smaller magnetic lever arm A. After some threshold, depending on the other MAES 
(turbulence) parameters, thermal effects become so important that super-A jets can be obtained 
with slow-rotators (a;^ < 1)- Right: Two new extreme MAES with an additional heating Q. If Q 
is very large at the disc surface (/ ~ 1), hot and dense jets are produced (solid curves) whereas 
cold and very tenuous jets (dashed curves) are obtained if Q is almost incxistcnt (/ = 5 10^^). 
Lower pannels show the effective turbulent heating Fe// and the prescribed entropy source Q 
at a radius Tq, normalized to the same quantity. 



This expression is consistent with global energy conservation, the parameter / being free 
and varying from ("cold" MAES) to 1 ("warm" or magneto-thermally driven jets). A 
vahic of / larger than unity would require an extra source of energy. We need now to 
specify the vertical (self-similar) profile of Q. Obviously, this introduces so many degrees 
of freedom that we do not dare anymore to look for the parameter space. Instead, we will 
look for extreme configurations and try to derive quantitative results. 



5.3.2 Dynamical effects 

The steady-state energy equation (79) can be explicitely written 

= prf = ,:5£_^.^^£,3^.Vr-S,.VP. (81) 

Dt Dt Dt 7 — 1 nip 

We then see that the main influence of a non-zero Q is on the temperature and pressure 
vertical gradients (thin disc), thus at the disc surface. But this is precisely the place 
where such gradients are required for ejection (Sect. 4.1.3)! Therefore, allowing some 
energy deposition at the disc surface has two major effects: 
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• (1) The initial jet temperature (T+) is increased and may thereby provide a non- 
neghgible initial enthalpy H ("warm" jets). Indeed, the Bernoulli integral becomes 



E{a) = — + H + - ^* 



rB^ 



V 




Q{s',a) 
,+ p{s' , a)up{s' , a) 



ds' 



(82) 



where s is a curvilinear coordinate along the magnetic surface and s"*" represents the 
disc surface. If Q remains positive above s+, it offers an additional energy reservoir 
for plasma. Moreover, the total energy feeding a magnetic surface, 



is affected by the heating that already occurred in the underlying layers. Such an 
heating may provide a ratio T^/To larger than unity, possibly allowing to relax the 
constraint on minimum field lines inclination. 

• (2) If Q is relevant in the upper resistive layers, then it will increase also the plasma 
pressure gradient. This will enhance the ejected mass flux (and lower the magnetic 
lever arm) and might therefore have tremendous consequences on jet dynamics. 

We used only one type of vertical profile Jq{s), changing the value of the parameter /. 
The chosen profile provides Q ~ inside the disc (a; < 0.5), an increase until a maximum 
value (fixed by /) around x < 1, then a decrease to zero (adiabatic behaviour) after 
roughly one scale height (see Casse & Ferreira 2000b). As / increases (all other MAES 
parameters being held constant), one goes from cold (T"*" < Tq), tenuous (/t < 0.02) jets 
from fast rotators {loa > 1) to hot (r+ Tq), dense (k > 0.1) jets from slow rotators 
{loa < 1)- Namely, the presence of some chromospheric heating allows a smooth transition 
from "cold" (purely magnetically-driven) to "warm" (magncto-thcrmaly driven) jets. 

Another class of "cold" solutions can also be designed. Indeed, if local cooling is 
not sufficient inside the disc {Q > for x < 1), the plasma pressure increases which 
provides the disc a stronger support against both tidal and magnetic compression. As 
a result, the magnetic field lines can be more bent than in the previous adiabatic or 
isothermal solutions. Such a large curvature hinders mass to be ejected {k may be as 
small as 10~^ and A ~ 400) but a vertical equilibrium can nevertheless be reached^^. 
Thus, entropy generation inside the disc removes the limits found on the "cold" parameter 
space described earlier. 

5.3.3 Global energy conservation 

We suppose that a MAES is settled around a central object of mass M, between an inner 
radius and some outer radius re- At this outer radius, the structure is fed with an 
accretion rate M^e = Ma{re)- Mass conservation then writes Mae — 2Mj = Mai and the 
fraction of ejected mass is 



the last expression being valid only for a very small ejection efficiency ^ <C 1. The local 
energy conservation equation is 




(83) 




(84) 




+ Smhd - u-T = pTup ■ VS - Teff 



(85) 



whereas the second law of thermodynamics (79) provides 



Q = pTUp ■ VS = Teff + {Tturb — ^turb) + ^ext — V • Srad 



(86) 
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Figure 6: Sketch of the volume used for the calculation of the global energy conservation. 



To get the global energy conservation, we integrate this equation on the volume V 
occupied by the disc. We thus define T,'^ and S~ as the disc surfaces^^ and Sj (Sg) the 
lateral surfaces at r = (r = re, see Fig. (6)). After integration, we get 

Pace "I" Pext "I" Pturb ~ "^Pjet "I" ^Prad (^'^) 

where the accretion power Pacc} i^* the power released by the accretion flow, is the 
difference between what comes in at re and goes out at rj. As said previously, we assumed 
Pext = Pturb = (neither an external source of energy, nor a significant input of turbulent 
energy). Thus, all available power Pace is shared between radiative losses at the disc 
surfaces P^-ad = J^± S^ad ' dS and jet power Pjet = /s± pUpE{a) ■ dS. It is usefull to 
introduce the fiducial quantity 

Pub = mib^^ ^2.510 ^—j [ ,,.rM^/yr ) I^TauJ ^''^ 

where 77/45 is a term roughly equal to unity. Within our self-similar framework, energy 
conservation of a thin (or slim) disc writes 



lib 



2Pjet ^ A 

Plib 1 + A 



5J 



qBo 



2 1 + A 



2Prad ^-^ j-^Pdiss Paec 2Pjet (91) 

Plib Plib Plib 

See Casse & Ferreira (2000b) for the derivation of these expressions. Three important 
remarks. First, the absolute limit for the ejection efficiency is ^ = 1; Second, the real 
energy release is Pace, comparable to Pub only for low ejection indices. Finally, the MHD 
Poynting flux feeding the jets depends directly on the amount of toroidal field at the 
disc surface. Thus, magnetically-dominated discs (A S> 1) may still produce some disc 
luminosity, provided \B'^\ < qBo (and / < 1). For / ~ we found solutions with a ratio 
Pjet/Prad varying between 0.1 and 10. 

5.3.4 Preliminary conclusions 

From this study, we are forced to conclude that thermal effects have an outrageous quan- 
titative importance on jet launching. In a way, this is fortunate since we can now recover 



^^Note that relativistic speeds are expected if such a MAES is settled around a compact object. 
^^The disc surface z = = x~^h{r) is precisely defined as the locus where Ur{r, z~^) = 0. 
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and understand the parameter space obtained with numerical MHD simulations. Al- 
though "cold", the simulated jets have enormous mass loads and correspondingly small 
magnetic lever arms ((eg. k ~ 0.6, A ~ 3.5, Ouyed &: Pudritz 1999). This may be easily 
achieved from an accretion disc where a significant entropy generation took place at the 
disc surface (/ > 0.1), followed right afterwards by a strong cooling (Q < for x > 1). 
In this way, very dense but cold (well, after some time) jets can indeed be produced. 
Remember that the bottom of the computational box has nothing to do with the real disc 
surface. 

A major drawback of such a conclusion is that any quantitative prediction (eg. the 
value of ^) requires to treat the energy equation. We must therefore inject our knowledge 
on MHD turbulence and its (possibly non local) energy transport properties ! Besides, 
illumination effects by the central protostar may also be important (since it heats up the 
disc surface). 



5.4 Observational predictions 
5.4.1 Accretion discs 

Whether or not the disc remains geometrically thin depends on its internal temperature, 
hence on the energy equation. Although the energy transport processes are unknown, 
an opaque disc must radiate at its surfaces all deposited energy. This translates into an 
effective temperature at the photosphere z = hphotif) 

aT^ff = V • Sraddz = (1 - /) [vm4 + v'mJp + Vv \rVn\^\ dz (92) 

The rhs of this equation can be estimated whereas the effective and central temperatures 
are crudely related through T,,f f ~ T^t^'^/'^, where r ~ K'Pohphot is the disc optical depth. 
Now, assuming that hph(,t{r) ~ h{r) and choosing the dominant disc opacity to be the 
grain opacity k = 0.1 T^^^ cm^g~^ (Bell & Lin 1994), allow us to estimate both the local 
disc aspect ratio e and optical depth r without even solving the energy equation. Once we 
have e = h/r, we can express all disc quantities in terms of the remaining free parameters. 
Hence, steady-state theory of MAES requires the following physical conditions in a disc 
settled around a low-mass protostar: 

^ = 5 10- «-V9(l - /)V9 f ^ V'V^^ 
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where Mgff = Mi(''o)/(l+-^)- Parameters /it ~ 1 and ^ are provided by the MHD solution. 
Note that the self-similar scaling h{r) oc r is here slightly violated since Ma oc r|. But 
such a deviation remains minor for ^ < 0.1. Care must be taken when looking at those 
expressions because A is imposed by Eq. (65) and depends thus on ie,am,'Pm,Xm)- But 
this choice allows to recover both radial scalings and values of a standard accretion disc 
when ^ = / = A = 0. To get a super-A solution with = am'PmfJ'^^'^ ^ 1 one probably 
needs to play around with / (for / = the parameter space is already known and requires 
large values of a^). 

Several observational diagnostics could reveal the presence of a magnetized disc driving 
jets: (i) the measure of a large (organized) disc magnetic field; (ii) optically thin regions 
or some lack of disc emission (for models with A ^ 1); (iii) a different spectral energy 
distribution (the effective temperature scales as riy' ^■*^^) for discs with a large ejection 
index ^. But usually, SED are tricky to interpret and we don't have the resolution yet to 
measure disc magnetic fields. 

5.4.2 Self-coUimated jets 

The most acurate tool to discriminate between models is to confront theoretical predic- 
tions with recent spatially resolved observations of the inner wind structure of TTS in 
forbidden emission lines of [O i], [S ii] and [N il]. Indeed, being optically thin, these lines 
carry information on both dynamic and thermodynamic properties of the whole volume 
of emission. So, one way to use this information is to construct the following synthetic 
observations and comparing them to real ones: 

• Emission maps, which must then be convolved to typical resolutions in order to 
predict what would be the observed jet morphology and collimation properties (eg. 
displacement of emission peaks, jet FWHM as a function of distance). 

• Line profiles, like those obtained using long-slit spectroscopy, and integrated pro- 
files. 

• Integrated line fluxes as well as their correlations with the disc accretion rate. 

• Forbidden line ratios, which reflect the values of the electron density, ionisation 

fraction and local temperature. 

While the first two observations offer constraints on the jet dynamics, the last ones test 
mainly the heating and ionization mechanisms along the jet. We already have a dynamical 
model providing us with jet density and velocity. The gas emissivity then requires to know 
the jet thermal and ionization states. With that in hand, we can easily compute synthetic 
observations and compare them to real ones. This has first been done by Safier (1993) 
using Blandford Sz Payne jet solutions. We can now do the same with models of MAES 
that take self-consistently into account the disc-jet interrelations. 

Cabrit et al. (1999), using parameterized temperature and ionization fraction, pro- 
duced synthetic maps and long-slit spectra that were nicely compatible with observations 
(Shang et al. (1998) did the same for X-winds). But more reliable predictions require 
actually solving for the jet thermal and ionization state, given some local heating mecha- 
nism. Our models being stationary, we cannot invoke shocks as being the heating process. 
Moreover, this would require the introduction of additional free parameters ! This is even 
worse if we rely on some small scale turbulence, providing a means to dissipate locally 
jet kinetic energy. The only self-consistent mechanism, intrinsic to MHD flows of low 
ionization, is ambipolar diffusion (see Sect. 2.2.3). 

Using this (unavoidable) effect to heat up jets, Garcia et al. (2001) solved the energy 
equation along the jet as well as its ionization state, taking into account several heavy 
elements (C, N, O, S, Ca, Mg, Fe...), photoionization heating and cooling by Hydrogen 
recombination lines. They used cold (isothermal) MAES solutions obtained by Ferreira 
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(1997) with ^ ~ 0.01. Jet widths and variations in hne profiles with distance and hne 
tracer are well reproduced. However, predicted maximum velocities are too high, total 
densities too low and, as a result, the low-velocity [O l] component is too weak. Thus, 
denser and slower MHD jets are required. 

In the solutions used, ejected plasma reaches its asymptotic speed ~ Cloro^/2X^^, 
which is typically 10 times the Keplcrian speed at the field line footpoint. Even though 
there is some inclination effect that decreases the observed jet velocity, it is still too large: 
emission from the jet inner region is important. Thus, one needs to decrease the jet 
terminal speed, which requires to diminish the magnetic lever arm A, ie. to increase the 
ejection index ^ (which automatically provides a denser jet). However, models with values 
of larger than 0.01 require an additional heating at the disc surface. How much exactly 
remains to be worked out. 

6 Some perspectives 

At this stage, one can have confidence in several things: (1) we know exactly, in a model 
independent way, how accretion discs can steadily drive jets; (2) jet properties (velocity, 
collimation) are strongly dependent on the mass flux allowed by the disc; (3) STich a mass 
flux is highly sensitive to the critical disc energy equation; (4) we don't know if real 
accretion discs ever meet the required physical conditions for MAES. 

The MAES paradigm is based on the belief that a large scale magnetic field is threading 
the disc on a large radial extension. This ensures some kind of cylindrical geometry 
(actually a conical Alfven surface), whereas ejection from a small disc region, as in X-winds 
(Shu et al. 1994), would provide a spherical expansion of the field lines (and a more or 
less spherical Alfven surface). Such a basic fact provides very different jet behaviours (see 
appendix B in Casse & Ferreira 2000a), but the underlying accretion/ejection mechanism 
remains exactly the same. As we saw, comparing synthetic to real observations is quite a 
powerful tool to eliminate MAES models and will probably help to discriminate between 
"disc winds" and "X-winds". We know now that, inside the MAES paradigm, YSOs 
require active chromospheres and coronae^^. Such a work must therefore be continued, in 
the hope that a characterization of their properties can indeed be achieved. In parallel, 
a thorough theoretical (analytical and numerical) work must be performed to understand 
how instabilities in magnetized discs may lead to turbulence and anomalous transport. 

What is not addressed by current models of MAES, is the interaction with the central 
object. However, including such an interaction and even allowing for some energetic events 
to occur there (eg. such as unsteady mass ejection) would not perturb the above picture. 
The main strength of the simplified version studied here (pure "disc winds" ) is to provide 
clear answers to the basic phenomena of correlated accretion and ejection. Anyway, 
such an interaction is unavoidable and, for a young protostar, may be the ingredient 
solving the "angular momentum problem". This problem arises from the idea that a 
protostar of masse M and radius R arises from the collapse of a rotating molecular cloud. 
The specific angular momentum of these clouds is typically 10^^ cm^/s (Goodman et 
al. 1993). A typical T-Tauri star (optically revealed protostar) rotates with a period 
of 8 days, which provides a specific angular momentum of the order of some 10^^ cm^/s. 
Where has all this angular momentum gone ? Evidently, we can assume that a (magnetic) 
braking occurred during the collapse, but no model provided yet a good account of the 
angular momentum transport. A conservative approach is to assume that, once formed, 
the protostar is rotating close to break-up, with a period P = 2t: yj E? /GM (~ 1 day). 
But the corresponding specific angular momentum is still an order of magnitude larger 

^^Kwan (1997) reached the same conclusion from energetic requirements on the low-velocity component of 
emission lines. 
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than that of a T-Tauri star. Thus, another very efficient braking must occur during the 
embedded phase. Moreover, there are some observational evidences that accretion discs 
play a role in braking down T-Tauri stars. Indeed, even though they are contracting and 
accreting material, T-Tauri stars seem to maintain a constant angular velocity (Bouvier et 
al. 1997). The idea is therefore to imagine an interaction with the accretion disc providing 
a negative torque on the protostar. 

If nothing prevents the accreting disc material to approach the stellar surface, an 
equatorial boundary layer will form at the interface (eg. Popham et al. 1993). But 
because matter is rotating faster than the protostar, such an interaction can only lead to 
a spin-up of the protostar. The cTirrcnt paradigm is a magnetospheric interaction with 
the disc. The protostar is believed to have a large scale magnetic field able to truncate 
the disc by forcing the accreting material to leave the equatorial plane and follow the 
magnetosphere. A spin-down is then obtained whenever the field is large enough so 
that the truncation radius (magnetopause) is located beyond the co-rotation radius Vco, 
defined by ^Igtar = ^disc(^co)- In this case, the protostellar angular momentum is deposited 
into the inner disc region... which then must expell it ! 

Some models rely on a very efficient disc turbulent viscosity that would radially trans- 
port both disc and protostellar angular momentum (eg. Li 1996). Or, if not throughout 
all the disc, just on a small region until it reaches a region located beyond the closed 
magnetosphere. There, open field lines could transport this excess of angular momentum 
in a jet ("X- winds" by Shu et al. 1994). As said previously, these X- winds carry away 
the exact angular momentum at a rate allowing accretion: they are therefore disc-winds, 
but produced from a tiny disc region (see also Fendt & Elstner 2000). It may be more 
economic to use the protostellar angular momentum and rotational energy to directly 
produce another jet component. This can be done whenever the protostellar magnetic 
moment has the same polarity as the disc magnetic field. A magnetic neutral line forms 
at the equatorial plane, where rcconnection of the two fields takes place. Above this 
reconnection site, a fraction of the accreting matter can be loaded onto newly opened 
protostellar field lines (Ferreira et al. 2000). Such Reconnection X- winds arise from a dif- 
ferent ejection mechanism: they expell disc material thanks to the protostellar rotational 
energy. This is the reason why they can brake down very efficiently the protostar, on time 
scales compatible with observations. 

It took almost 30 years to prove that accretion discs can indeed launch jets that carry 
away a sizeable fraction of their angular momentum. Investigations of the magnetospheric 
interaction between a protostar and its circumstcllar disc are just beginning and the 
situation is far more complex than for discs. Lots of analytical efforts have still to be 
done, probably guided by the insights provided by numerical simulations. 
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